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Abstract 

A small electromagnetic sampling calorimeter, installed in the ZEUS experiment 
in 1995, significantly enhanced the acceptance for very low x and low Q 2 inelastic 
neutral current scattering, e + p — ► e + X, at HERA. A measurement of the proton 
structure function Fi and the total virtual photon-proton (7*2?) cross-section is 
presented for 0.11 < Q 2 < 0.65 GeV 2 and 2 x 10~ 6 <i<6x 10 -5 , corresponding 
to a range in the 7*2? cm. energy of 100 < W < 230 GeV. Comparisons with various 
models are also presented. 
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1 Introduction 



An early discovery at HERA was the rapid rise of the proton structure function, F 2 (x, Q 2 ), 
as the Bjorken scaling variable, x, decreases at low x for photon virtualities Q 2 > 10 GeV 2 
[U 0. The ZEUS [[| and HI Collaborations have extended the measurement of F 2 
down to a Q 2 value of 1.5 GeV 2 . One of the most interesting features of the recent 
data is the persistence to the lowest Q 2 of the rapid rise of F 2 with decreasing x. The 
predictions of Gliick, Reya and Vogt (GRV) which result from the dynamic generation 
of parton densities via next-to-leading order (NLO) perturbative QCD (pQCD) DGLAP|| 
evolution of valence type distributions starting at a very low scale, Ql = 0.34 GeV 2 , are 
in broad agreement with this observation. It is surprising that leading twist NLO pQCD 
can describe the data down to Q 2 = 1.5 GeV 2 . Other global fit analyses based on NLO 
pQCD, such as those performed by MRS |7j and CTEQ ||, typically have much higher 
starting scales Qq = 3 - 5 GeV 2 . It then becomes an interesting question to determine at 
which Q 2 the behaviour of F 2 becomes dominated by non-perturbative contributions. 

F 2 is related to the total virtual photon-proton (7*2?) cross-section by o~loi ~ {^ 2 ot/Q 2 )F 2 
for x <C 1. At fixed Q 2 , the rapid rise of F 2 with decreasing x is equivalent to a rapid rise 
of the total j*p cross-section with cm. energy, W (W 2 ~ Q 2 /x in this kinematic regime). 
At high W, oJot can be described by a power law behaviour, aJot °^ W 2X . For Q 2 > 1.5 
GeV 2 , the power A is between 0.15 and 0.35 || £§]. In contrast, the total cross-section for 
real photon-proton scattering (photoproduction, with Q 2 = 0) shows only a modest rise 
with W, A = 0.08 @, consistent with the energy behaviour of the total pp cross-section. 
Regge theory has been used successfully by, e.g., Donnachie and Landshoff (DL) || to 
describe the energy dependence of the total cross-section for hadron-hadron scattering 
and real photon-proton scattering, but their prediction (|Tt| for virtual photon-proton 
scattering fails to describe the data for Q 2 > 1.5 GeV 2 . Different groups (CKMT [|TT|1 , 
BK[fr2[, ABY|T3|], and ALLM[]rjJ) have used a variety of approaches to connect the very 
low Q 2 behaviour with high Q 2 pQCD. CKMT extend the Regge prediction by including 
Q 2 dependent absorptive corrections that modify the effective pomeron intercept, resulting 
in a Q 2 dependent A up to some Ql in the range 1 to 5 GeV 2 ; pQCD is then used to 
evolve from this to higher Q 2 . BK describe the Q 2 behaviour using a generalised vector 
dominance model (GVDM): the low Q 2 region is controlled by the contributions of the low 
mass vector mesons, and the higher mass contributions are adjusted to provide agreement 
with pQCD predictions using a standard set of structure function parametrisations at 
larger Q 2 . ABY extend their high Q 2 QCD-inspired parametrisation into the low Q 2 
regime, and modify the evolution of a s so that it saturates at a finite value. ALLM 
introduce parametrisations that interpolate between the Regge and pQCD regimes. A 



review is given in reference [[15 



To study the transition from the hadronic type behaviour at Q 2 ~ to the deep inelastic 
scattering (DIS) regime (Q 2 ^> 1 GeV 2 ), the kinematic coverage of the ZEUS detector was 
substantially extended starting in 1995 with the installation of a new beampipe calorimeter 
(BPC)Q Here we report on the measurement of F 2 and a] ot p for 0.11 < Q 2 < 0.65 GeV 2 
from e + p scattering at a/s = 300 GeV using the ZEUS detector with the new BPC. This 
analysis is based on 1.65 pb _1 of data taken during the 1995 HERA run. Recently the HI 



1 The new BPC replaced a previous calorimeter described in [n6| 
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Collaboration has reported an F2 measurement in four bins in the Q 2 range of 0.35 to 0.65 
GeV 2 in a somewhat different W range [Jl7| . For Q 2 > 0.23 GeV 2 , the E665 Collaboration 
reported a measurement of the proton structure function at much higher than this 

analysis. 



2 Kinematics 

Inelastic positron-proton scattering, 

e + p e + X (1) 

can be described in terms of two kinematic variables, x and Q 2 , where x is the Bjorken 
scaling variable and Q 2 the negative of the square of the four-momentum transfer. In the 
absence of initial and final state radiation, Q 2 = —q 2 = —(k — k') 2 and x = Q 2 /(2P ■ q), 
where k and P are the four-momenta of the incoming positron and proton respectively 
and k' is the four-momentum of the scattered positron. The fractional energy transfer 
to the proton in its rest frame, y, can be related to x and Q 2 by y = Q 2 /(sx), where 
s is the square of the e + p cm. energy which is given by s ~ 4E e E p . Here, E e (27.5 
GeV) and E p (820 GeV) are the positron and proton beam energies, respectively. The 
kinematic variables, y and Q 2 , are related to the energy, E' e , and angle with respect to the 
proton beam direction, 9 e , of the scattered positron. We also use $ = ir — 9 e , the angle 
with respect to the positron beam direction. Scattering at low Q 2 results in positrons 
emerging at small 

y = 1 - Jkl + costf) « 1 - §, Q 2 = 2E e E' e {l - costf) « E e E'J} 2 . (2) 



3 Experimental setup and kinematic reconstruction 



The ZEUS detector |19[] is a general purpose magnetic detector at the HERA collider. To 
enhance the acceptance of the detector at small Q 2 , two beam pipe calorimeter modules 
(BPC) |21J located on two sides of the beam pipe at 2937 mm from the interaction point 
in the rear (positron) directionP] were installed, as shown in figure |l](a). The BPC covers 
positron scattering angles relative to the incident direction of 15 to 34 mrad. At these 
small angles, the maximum possible scattered positron energy is equal to the beam energy, 
27.5 GeV. The beam pipe has two low-mass windows (0.016 radiation length (r.l.)) in 
front of the BPC to allow positrons to exit the beam pipe with minimal interference. 

The BPC is an electromagnetic scintillator sampling calorimeter. The passive absorber 
layers consist of twenty-six 13.8 cm x 13 cm x 0.92 r.l. thick tungsten alloy plates; the 
active layers consist of 7.9 mm wide and 2.6 mm thick scintillator strips alternating each 



2 The ZEUS right-handed coordinate system is defined with the origin at the nominal interaction point 
by the Z axis pointing in the proton beam direction and the X axis horizontally pointing towards the 
center of HERA. 
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layer in the horizontal and vertical directions. The scintillator strips are read out from 
one end using wavelength shifters (WLS). Each WLS is coupled to a miniature photo- 
multiplier tube (Hamamatsu R5600-03). The vertically oriented scintillator strips provide 
the X position measurement and the horizontally oriented strips the Y measurement. 
The readout electronics are similar to those used for the main ZEUS uranium scintillator 
calorimeter (CAL) [|T^]. The alignment is known to an accuracy of 0.5 mm from an 
optical geometrical survey. The distance between the two calorimeter modules on either 
side of the beam pipe is mechanically constrained to within 0.2 mm. Due to synchrotron 
radiation from the positron beam, the modules are placed asymmetrically around the 
beam. Only one of the two BPC modules is used for physics analysis due to the very 
small acceptance of the other module, which is used exclusively for alignment purposes. 
The typical geometrical acceptance is 10%. 



3.1 Detector simulation and response 



The BPC simulation is based on the GEANT [21] program, with an independent check 



performed using EGS4 [22]. The energy spectra for 1 to 6 GeV incident electrons in the 
simulation are in good agreement with test beam data taken at these energies, and are 
consistent with the design energy resolution of \l%j\fE. The non-linearity is found to 
be less than 1% in the simulation for 2 to 6 GeV incident electrons, in agreement with 
test beam data. 

The energy calibration was performed in situ using kinematic peak (KP) events^ in two 
steps: a relative strip-to-strip calibration, followed by an overall energy scale calibration. 
Figure |l](b) shows the fractional deviation of the KP energy from the mean value as 
a function of the scattered positron X impact position at the BPC after the relative 
calibration. The energy scale is uniform to within 0.5% over the BPC fiducial region, which 
extends to 8 mm from the edge of the BPC. The overall energy scale was established by 
comparing the KP energy spectra of data and simulation, which included QED radiative 
corrections. A x 2 was calculated between the two spectra, and the energy scale of the 
data adjusted to minimise it. 

The absolute energy scale obtained with KP events was checked using elastic p° events, 
e + p — > e + p°p. The scattered positron energy and position were measured in the BPC, 
and the three-momenta of the two p° decay pions were measured using the ZEUS Cen- 
tral Tracking Detector (CTD)p4|. Using the four-momentum of the p°, the scattered 



positron energy can be independently determined. Figure 0(c) shows the ratio of the 
measured positron energy to that determined from the CTD, for both simulation and 
data. Radiative corrections are responsible for the tail at low values and the fact that 
the distributions peak below unity. The agreement between the data and the simulation 
is very good. Consequently, it was concluded that the energy scale determined using the 
KP events was accurate to within 0.5%, and that the BPC resolution was well modeled 
in the simulation. 

The BPC was located only 4.4 cm from the beam in the horizontal direction, and received 



3 A cut yjB < 0.04 (see below) selects scattered positrons whose energy distribution sharply peaks 
within 2% of the beam energy, providing a good calibration source B3|. 
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12 kGy of radiation during the 1995 HERA running period. The resulting damage caused 
a drop in the energy scale of up to 2.5% for the regions closest to the beam, determined 
using KP events. To correct for this, the data were separated into 4 time periods and the 
energy calibration performed separately for each period. The amount of radiation damage 
to the BPC was also determined using a movable cobalt source scan calibration system 
2H] and by measuring the response of a sample of scintillator strips after disassembling 



the calorimeter at the end of 1995 after data taking. The degradation of each individual 
scintillator strip was built into the EGS4 simulation, which showed that the non-linearity 
from 10 GeV to full scale due to radiation damage was less than 1%. 

The BPC measured the arrival time of the positron with an accuracy determined to be 
0.4 ns for positrons with an energy greater than 6 GeV. 



3.2 Positron identification and position reconstruction 

The scattered positron position in the BPC was determined using the logarithmic energy- 



weighted shower position |26| using scintillator strips containing more than 4% of the total 
shower energy. A RMS resolution of 1.3 mm with a maximum systematic shift of 0.5 mm 
was obtained for 5 GeV incident electrons in the beam test. This was well-reproduced 
in the simulation. The position resolution improves gradually with increasing positron 
energy. The resolution was determined from the simulation to be 0.6 mm at 27.5 GeV for 
positron impact positions within the fiducial region of the BPC. 

As a cross-check of the absolute position, QED Compton events e + p — > e +r yp were used. 
Since both the photon and the scattered positron are detected in the BPC modules, 
these events provide an over-constrained kinematic system. Using the accurately known 
distance between the two BPC modules and the QED Compton event kinematics, the po- 
sition of each BPC module relative to the positron beam was determined. The agreement 
with the survey is better than 0.5 mm. The accuracy of the QED Compton method is 
dominated by the uncertainty of the calorimeter energy scale: 0.5%, corresponding to 0.5 
mm. 

Positron identification was performed using the transverse size of the shower. The second 
moments of the shower in the X and Y directions, ax and cry, were calculated using the 



logarithmic energy weighted method mentioned above. The combination y \a x + er y )/2 
was required to be less than 7.5 mm. This yielded a positron acceptance in excess of 95% 
at 7 GeV, rising to 98% above 12 GeV, while rejecting hadrons and positrons that have 
preshowered in the beam-pipe wall, which have a much wider transverse width. Using a 
sample of KP positrons, the transverse energy profile in the simulation was tuned to that 
of the data. 



3.3 Vertex determination and luminosity measurement 

The position of the e + p interaction vertex is needed to determine the positron scattering 
angle. The Z position of the vertex was measured using the CTD on an event-by-event 
basis with a typical resolution of 4 mm. For events with no CTD vertex information 
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(about 8% of the total), the Z position of the vertex was set to the average position of 
the full data sample, <Z> = —3 cm. The longitudinal size of the luminous region was 
about 12 cm (r.m.s.). The mean values of the X and Y vertex positions, determined on 
a run-by-run basis, were used. The transverse sizes of the beam in X and Y were about 
300 /im and 70 jum, respectively. 

The luminosity was measured via the bremsstrahlung process e + p — > e + 7p, using a 
separate electromagnetic calorimeter detector system (LUMI) |[27|| . A lead-scintillator 
calorimeter positioned at Z = —107 m, accepting photons with scattering angles less than 
0.5 mrad, was used for the luminosity measurement. The uncertainty of the luminosity 
measurement for the data sample used in this analysis is 1.5%. A second electromagnetic 
calorimeter, positioned at Z = —35 m, was used for tagging positrons in background 
studies. 



3.4 Reconstruction of the kinematic variables 

In this analysis, the kinematic variables y and Q 2 are reconstructed using the energy and 
angle of the scattered positron, determined using the BPC and the CTD vertex position 
with equation Using this method of reconstruction ("electron method"), y and Q 2 can 
be determined reliably over the kinematic range y > 0.1 and Q 2 > 0.1 GeV 2 . The y 
resolution is 0.02 to 0.04, and the Q 2 resolution is 6 to 8%. 

As the ZEUS CAL is an almost hermetic detector, it can be used to measure the hadronic 
system, denoted by X in equation [[]. The following quantities are reconstructed using the 
CAL, 

$cal = ^Ei(l - cos6»j), y JB = Etot = ^E h (3) 

where Ei is the energy measured in the i th CAL cell and 0, is the polar angle of the center 
of the i th calorimeter cell with respect to the positive Z axis; the sum extends over all 
cells in the CAL. The quantity yjB provides a measure of the kinematic variable y, and 
has superior resolution at very low y compared to that from the electron method. E tot 
is the measured total energy of the hadronic system. These quantities were used in the 
trigger and event selection for rejecting background, reducing QED radiative corrections 
and controlling event migration at low y. 



4 Trigger, data taking and event selection 



ZEUS selects events online using a three-level trigger system ||28|| . Both the energy and 
timing information from the BPC were used for the First Level Trigger (FLT). An energy 
cut of 6 GeV was made and the timing was required to be consistent with an e + p collision. 
Proton-gas events occurring upstream were also rejected by timing measurements made by 
scintillation counter arrays situated along the beamline at Z = —730, —315, and —150 cm, 
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respectively. For the Second Level Trigger, CAL timing information was used to reject non 
e + p events. An approximate value of yjB, determined from the CAL energies assuming 
an interaction vertex at Z = 0, was required to be greater than 0.02 and the total CAL 
energy was required to be greater than 3 GeV. No additional BPC cuts were imposed at 
the Third Level Trigger. 

The FLT efficiency was studied using a sample of events triggered only by the CAL. All 
of the offline event selection cuts (see below) were applied to this sample. The trigger was 
found to be fully efficient for BPC energies greater than 7 GeV, as shown in figure |l|(d). 

The offline event selection cuts are as follows. The reconstructed positron is required to 
be within the BPC fiducial region, to have more than 7 GeV energy, and to pass the 
shower width cut described in section |3.2j . The BPC time is required to be within 3 ns of 
the time for e + p interactions. If the event vertex is well reconstructed with the CTD, the 
Z position of the vertex is required to be within the window -40 cm < Zyrx < 100 cm. 
The quantity 5 = 5c al + 2E' e , is required to lie in the range 35 to 60 GeV; 6 equals twice 
the positron beam energy (55 GeV) for a completely contained event, but the distribution 
peaks at much lower values for photoproduction events where the scattered positron is 
lost in the rear direction. A cut yjB >0.06 reduces migration of events from very low y, 
where the resolution of the electron method is poor. Finally, if timing information from 
the CAL is available, the time is required to be consistent with an e + p collision. After 
cuts, 109105 events remain in the data sample. 



5 Analysis 

5.1 Physics simulation 

A physics simulation is used to determine the radiative corrections and the acceptance 
of the detector. The starting point for the simulation of e + p collisions in the y and Q 2 
region of this measurement is the program DJANGO 6.22, which interfaces the programs 



HERACLES [[29| 4.5.1, ARIADNE gg 4.06 and LEPTO [|3TJ 6.4.1. The HERACLES 
program calculates the structure functions F 2 and Fl from an input set of parton density 
functions. From these, it calculates the differential cross-section including initial and 
final state radiation, and the full one-loop virtual corrections. ARIADNE implements the 
colour dipole model for gluon radiation between the struck quark and the proton remnant. 



Finally, LEPTO handles the fragmentation using the program JETSET [32] . 



Several modifications were made to the program to generate events at low Q 2 . (1) The 
Donnachie-Landshoff parameterisation [1C] was used to calculate F 2 in the Q 2 range of this 



measurement . 

(2) The longitudinal structure function Fl was set equal to zero at this stage. (3) Elastic 
vector meson (VM) p° events were generated. The VM events were re-weighted accord- 
ing to the cross-section and W dependence recently measured at ZEUS using the BPC. 
This contribution amounted to 6% of the events ||33jl - (4) Diffractive events were gener- 
ated according to d 2 a/dtdM x ~ e bt /(M x + Q 2 - M 2 ) 1 - 1 0, where t is the square of 
the four-momentum transferred to the outgoing proton, Mx is the invariant mass of the 
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hadronic final state, and 6 = 6 GeV~ 2 . The fraction of diffractive events (around 25%) 
was determined from the data by counting the events with a characteristic rapidity gap; 
that is, a region of little or no hadronic activity between the forward edge of the CAL 
and the jet from the struck quark. The overall acceptance is weakly sensitive to changes 
in the relative contributions of the different event types (see section 5.5). 

A generated event sample in the region Q 2 > 0.05 GeV 2 and y > 0.03 corresponding 
to more than twice the luminosity of the data was passed through the complete ZEUS 
simulation chain, which is based on the GEANT |H] program, and then processed using 
the same offline reconstruction software as for the data. The general characteristics of the 
data are well-described by the simulation, as shown in figure ^|(a-c). 



5.2 Binning of the data 

The data are binned in the variables y and Q 2 , which makes efficient use of the available 
phase space, as the lower electron energy cut corresponds to an upper y cut. The accessible 
kinematic region lies between y values of 0.08 and 0.74, and Q 2 values of 0.1 and 0.74 
GeV 2 . The sizes of the bins are chosen based on the experimental resolution of the 
kinematic variables and the number of events. At low y, the bin widths in y are chosen to 
be twice the y resolution; for y greater than 0.37, bins of approximately constant width 
are used. The lowest Q 2 bin has a width 2.5 times the Q 2 resolution; higher Q 2 bins have 
a constant width in log(Q 2 ) which results in approximately constant numbers of events 
in each Q 2 interval. 

For positrons within the BPC fiducial region, the efficiency of the event selection cuts is 
close to 100% for y < 0.5, decreasing to about 70% at y = 0.7. The purity, defined as the 
fraction of events reconstructed in a bin that were generated in that bin, is typically 50%. 



5.3 Background determination 

The background from beam-gas interactions is determined using unpaired positron and 
proton bunches. The size of this background is 1%, and is subtracted statistically. The 
dominant source of background comes from photoproduction interactions where the scat- 
tered positron escapes through the rear beam pipe and a fake positron is reconstructed 
in the BPC. Typically such events have much lower 5 values than signal events. Photo- 
production background events were generated using the PYTHIA program |35| with the 



ALLM cross-section parameterisation [I4J. This PYTHIA sample is used to perform a 
bin by bin subtraction of the photoproduction background. Figure ||(d) shows the 5 spec- 
tra for data (solid circles), signal simulation (dashed line), photoproduction simulation 
(shaded region), and the sum of the signal and photoproduction simulations (solid line). 
The good agreement between this sum and the data suggests that the photoproduction 
background is well simulated, and gives a contamination of a few percent in most bins, 
rising up to 15% in the highest y bins. 

As a cross-check, use was made of the very small angle LUMI positron detector to measure 
the photoproduction background directly. This detector accepts scattered positrons with 
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Q 2 < 0.01 GeV 2 and 0.2 < y < 0.6 and may be used to tag events with a fake BPC 
positron signal that pass the event selection cuts. The measured background is shown as 
the triangular points in figure @(d); the y < 0.6 cut limits the measurement to 5 < 35 GeV. 
Once again, the measured points are in good agreement with the simulated background. 

The contamination from higher Q 2 DIS events is less that 0.1% based on searches in the 
CAL for additional positron candidates in both data and the simulation. 



5.4 Determination of F 2 and cr] f 

The double-differential e + p cross-section for inelastic scattering can be expressed in terms 
of the total cross-section for virtual transverse (T) and longitudinal (L) photons: 



d 2 a(ep — > eX) 
dydQ 2 



a T (y, Q 2 ) + ea L (y, Q 2 ) (1 + S r (y, Q 2 )) (4) 



where the flux V « a(l + (1 — y) 2 ) / (2irQ 2 y), the fractional flux of longitudinally polarised 
virtual photons e w 2(1 — y)/(l + (1 — y) 2 ), and 5 r is the radiative correction factor. 
These cross-sections can be related to the proton structure functions F 2 and F L by F 2 = 
(Q 2 + Am 2 p x 2 ){a T + a L )/(A'K 2 a{l-x)), F L = (Q 2 + 4m 2 x 2 )a L /(4:ir 2 a(l — x)) and the total 
virtual photon-proton cross section by oJJ: = cr T + <j L . In the Q 2 range of this analysis, 
the contribution from Z° exchange is negligible. The radiative correction to the Born 
cross-section, S r , is a function of y and Q 2 , but to a good approximation is independent 
ofF 2 . 

An iterative procedure is adopted to extract the sum o T + eC r L . Monte Carlo events are 
generated as described in section |5.1| to determine the acceptance and efficiency of the 



cuts in each bin. The first ax + eo"L values are then fit with a simple functional form, 
inspired by the ALLM parameterisation ||14|| . The result of the fit is used to reweight the 



input distributions in the simulation event by event to re-evaluate the acceptance and 
efficiency New values of a? + ecr^ in each bin are calculated and the procedure repeated 
until the change between subsequent iterations is less than 0.5% in all bins; three iterations 
are required. The relative fractions of diffractive and VM events are kept fixed in this 
procedure. 

Once the ax + ^cr L values are determined, F 2 and a]^ are calculated assuming o L to be ei- 
ther zero or the value given by the Vector Dominance Model (VDM), <tl = K(Q 2 /My)<jT- 
My is set equal to the mass of the p° (0.77 GeV) and K to 0.5. 



5.5 Systematic uncertainties 

To estimate the systematic uncertainties of the measured F 2 values, the following checks 
were performed. In each case, an aspect of the event selection, reconstruction of kinematic 
variables, or F 2 determination was modified, the procedure described above was repeated 
and the change in F 2 noted. (1) The 0.5% uncertainty on the energy scale (see Section |3.1| ) 
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gave an effect of approximately 3% for F 2 ; the non-linearity of the BPC, estimated to be 
less than 1%, resulted in changes of less than 5%. (2) Varying the absolute position by 
0.5 mm produced changes of less than 6%. (3) To estimate the uncertainty due to the 
electron finding efficiency the shower width cut was varied by 1 mm yielding a change up 
to 2% at high y. (4) The uncertainty due to the CAL event selection cuts was checked by 
varying the CAL energy scale by 3%, changing the CAL minimum cell energy threshold, 
varying the 5 cut by 2 GeV and varying the yjB cut by 0.01. The effect on F 2 is small 
for moderate y bins, and reaches 6% at high and low y. (5) The estimated uncertainty 
on the photoproduction background was 30%, obtained from a comparison of the various 
methods described above. Consequently, the amount of background subtracted was varied 
by the same amount. The effect on F 2 increases with y, reaching 4% in the highest y bins. 
(6) The fraction of diffractive and VM events was varied by 25%. The effect on F 2 was 
small except for the lowest y points where it reached 4%. (7) Uncertainties due to the 
description of the hadronic final state in the simulation were estimated by comparing the 
results from different simulation programs(EPSOFT PYTHIA |J§, HERWIG 



constraints were provided by comparing hadronic distributions measured in the data with 
those predicted by the simulations. An error of 2% in F 2 for y less than 0.4, rising linearly 
to 5% at y = 0.74, was assigned. (8) The uncertainties in the radiative corrections, which 
modify the Born cross-section by 10-15%, were investigated. As mentioned in section 
|3TT| , the tail of the Ebpc / E ca i c distribution for elastic p° events (figure 1(c)) is due to 
initial state radiation. Comparison of the simulation with data in that region provides 
an estimate of the uncertainty in the radiative correction. This uncertainty represents 
possible changes of 3 to 4% on F 2 . 

The extracted uncertainties after each of the above systematic checks are displayed in 
table |A| (see Appendix). The total systematic error for each bin was determined by 
adding the changes to F 2 from different checks in quadrature. These are shown in table [l]. 
The systematic errors are around 6% for moderate y bins and are dominated by BPC 
calibration and radiative correction uncertainties. For low y, uncertainties in the CAL 
energy scale, as well as the dependence of the acceptance on the fraction of diffractive 
events leads to uncertainties up to 11%. At high y, uncertainties on the contributions 
from photoproduction background and the properties of the hadronic final state resulted 
in errors as high as 8%. 

The uncertainties in the luminosity measurement and trigger efficiency contribute to the 
overall normalisation error, which amounts to 2.4%. The statistical errors are in the range 
2-4%. 



6 Results 

In the following, F 2 is presented as a function of x (or y) and Q 2 , and o~l i as a function 
of W 2 and Q 2 . While the former vanishes as Q 2 — > 0, the latter is expected to extrapolate 
smoothly to the total photoproduction cross-section. 

The F 2 results are displayed in table [I] in bins of Q 2 and y and plotted in figure |3] as a 
function of x for different Q 2 bins. Here Fl is assumed to be zero. Assuming Fl as given 
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Table 1: Measured F 2 values, with the assumption Fi = 0. The first error indicates the 
statistical error, and the second indicates the systematic error. The values in the column 
marked AF 2 [Fl\ show the change in F 2 assuming a value for Fl given by VDM (see text), 
and are not included in the systematic errors. The bin boundaries in y are 0.08, 0.16, 
0.23, 0.30, 0.37, 0.45, 0.54, 0.64, and 0.74 and in Q 2 0.1, 0.13, 0.17, 0.21, 0.27, 0.35, 0.45, 
0.58, and 0.74 GeV 2 . 



by VDM, the effect on F 2 is typically around 1-2% for most bins, and increases F 2 by up 
to 7% in the lowest x bins, as shown in table |1| 



Also shown in figure |3| are the data from the E665JTBJ experiment at similar Q 2 , but much 
larger x values, and four recent points from the Hl|17| experiment with Q 2 < 0.65 GeV 2 . 



One observes a rise of F 2 by a factor 1.5 to 2 from x near 10~ 3 to x around 10 -5 . In the 
bottom of figure |3|, the F 2 values for Q 2 =1.5, 3.0 and 6.5 GeV 2 are shown, taken from 
previous HI and ZEUS publications P, f|, [L7] an d from E665. They illustrate the rapid 
rise of F 2 with decreasing x as Q 2 increases. 

Curves from various theoretical models are overlaid. In general DL and CKMT predictions 
are 15% and 10% lower than the data respectively, while the BK prediction is 15% higher. 
No attempt was made to modify the parameters of these models to fit the data. The value 
of F 2 given by GRV is approximately 35% of the measured value at Q 2 = 0.4 GeV 2 , rising 
to about 80% at Q 2 = 0.65 GeV 2 . At larger Q 2 the GRV predictions reproduce the rapid 
rise of F 2 , but tend to lie somewhat above the data. The ABY parametrisation, which 
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included preliminary ZEUS BPC F 2 results in the parameter fit, gives a good description 
of the data. 

Figure |] shows old at low Q 2 and at higher Q 2 taken from previous ZEUS analyses and 
from the HI and E665 experiments as a function of W 2 . The total cross-section for real 
•yp scattering is also shown |3~9|] . The curves in figure f| show the expectations 

from the soft Pomeron model of DL (dotted curve) and the pQCD model of GRV (dashed 
curve). The W dependence of the DL predictions is given by old oc W 2X with A = 0.08 
independent of Q 2 . The GRV model predicts a stronger variation with Q 2 and W 2 . 
Although the DL model curve is below the ZEUS data, its slope is in broad agreement 
for Q 2 < 0.65 GeV 2 ; however, this is not true for Q 2 > 1.5 GeV 2 . For Q 2 > 1.5 GeV 2 , 
the GRV model produces a rapid rise of old with W 2 that is in better agreement with 
the data. 



7 Conclusions 

In previous studies, ZEUS and HI have shown that for Q 2 > 1.5 GeV 2 and x <C 10 -2 , the 
proton structure function F 2 rises rapidly as x decreases, in agreement with models based 
on perturbative QCD. In this paper we have presented a measurement of F 2 and old for 
0.11 <Q 2 < 0.65 GeV 2 and 2 x 10~ 6 < x < 6 x 10~ 5 (100 < W < 230 GeV), covering 
the region between deep inelastic scattering and photoproduction. Similar results have 
recently been reported by Hl|l7|. In combination with data from E665 for x > 10~ 3 , F 2 
exhibits a modest rise in this low Q 2 region. Together with previously published HERA 
data for Q 2 > 1.5 GeV 2 , where the rise is more rapid, our results suggest that pQCD 
calculations can account for a significant fraction of the cross-section starting at Q 2 ~ 1 
GeV 2 . 
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Figure 1: (a) Schematic layout of the BPC and the beampipe. The BPC is located at 
Z=-2937 mm, and the inner edge is at X=44 mm. (b) The fractional deviation from the 
mean energy for KP events as a function of the scattered positron X impact position at 
the BPC. (c) A comparison between MC and data of the ratio of the measured positron 
energy in the BPC, E B pc, to the calculated positron energy, E ca ; c for elastic p° events, 
(d) The FLT trigger efficiency as a function of BPC energy. 
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Figure 2: Comparisons between simulated and data distributions for the variables: (a) 
scattered positron energy; (b) Q 2 ; (c) x; and (d) 5 . Data are shown as solid circles, 
photoproduction simulation as shaded regions, and the sum of the signal and photopro- 
duction simulations as solid lines. In figure (d) the signal simulation is shown as a dashed 
line and the measured background as triangular points. The structure function in the 
simulation has been reweighted to that measured in the present analysis. 
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Figure 3: F 2 (x, Q 2 ) as a function of x for different Q 2 bins. The data from this analysis, 
ZEUS BPC 95, are shown as solid dots, with E665, HI and previous ZEUS points shown 
as open squares, open circles and solid triangles, respectively. New points from H1[|I7| at 
low Q 2 are also shown as solid inverted triangles (the point at Q 2 = 0.35 GeV 2 is displayed 
in the 0.4 GeV 2 panel; additional points at 2.5 and 3.5 GeV 2 not shown.) Curves from 
the models of DL, CKMT, BK, ABY and GRV are overlaid. 
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Figure 4: The total virtual photon-proton cross-section alJ as a function of Vr 2 . The 
data from this analysis (ZEUS BPC 95), previous 1994 ZEUS analyses, HI, and E665 
are shown. The total cross-section for real photon-proton scattering from ZEUS, HI and 
photoproduction experiments at low W are also shown. The predictions of DL and GRV 
(at the ZEUS Q 2 values) are indicated by the dotted and dashed curves respectively. 
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Appendix A: individual effects of systematic uncertainties 
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Table A: The determination of the systematic errors was discussed in section 5.5. Q 2 , y, and measured F 2 values are listed in columns 
1 to 3. The individual effects of each of the systematic checks on F 2 are displayed (in percentages of F 2 ) in columns 4 to 27. 



